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Abstract 

We propose an explicit way to generate a large class of Operator scaling Gaus- 
sian random fields (OSGRF). Such fields are anisotropic generalizations of self- 
similar fields. More specifically, we are able to construct any Gaussian field 
belonging to this class with given Hurst index and exponent. Our construction 
provides - for simulations of texture as well as for detection of anisotropics in 
an image - a large class of models with controlled anisotropic geometries and 
structures. 
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1. Introduction 



Random fields are a useful tool for modelling spatial phenomena such as 
environmental fields, including for example, hydrology, geology, oceanography 
and medical images. Particularly important is the fact that in many cases these 
random fields have an anisotropic nature in the sense that they have different 
geometric characteristics along different directions (see, for example, Davies and 
Hall (0), Bonami and Estrade ([4]) and Benson, et al.(3)). 

Moreover, many times the model chosen has to include some statistical de- 
pendence structure that might be present across the scales. For this purpose, 
the usual assumption of self-similarity is formulated. Unfortunately, the clas- 
sical notion of self-similarity (see [15|), defined for a field {X{x)}^^g^d on R'^ 

by 

for some H E M. (called the Hurst index), is genuinely isotropic and therefore 
has to be changed to fit anisotropic situations. 



For this reason, there has been an increasing interest in defining a suitable 
concept for anisotropic self-similarity. Many authors have developed techniques 
to handle anisotropy in the scaling. The main papers that have to be men- 
tioned in this context are those of Hudson and Mason, Schertzer and Lovejoy 
(see 



]J,ll7|,|l8|). 



This motivated the introduction by Bierme, Meerschaert and Scheffler of 
operator scaling random fields (OSRF) in Q. These fields satisfy the following 
scaling property : 

{X(a^x)},e„. = {a"X{x)}^^^, , (1.1) 
for some d x d matrix E with positive real parts of the eigenvalues. 

A large class of random fields obeys this property. For example the Fractional 
Brownian Field (FBM) and the Fractional Brownian Sheet (FBS) are both Oper- 
ator Scaling Gaussian Random Fields (OSGRF) with exponent E = Id. Denote 
< •, • > the Euclidean scalar product of defined for any x = {xi, ■ ■ • , Xd) G K'' 
and y — (yi, ■ ■ ■ ,yd) £ M.'^ a.s < x,y >= J2i=i ^iVi- Recall that the FBS is the 
Gaussian field [Bh^, -- ,HA^)}xm'^ defined for some {Hi, • • • , Hd) G (0, 1)'' as : 

(e'<^'«> - l)dW{i) 



where dW is the Fourier transform of white noise on WK This Gaussian field 
enjoys with the following scaling property : for all (ai, ■ • • , ad) G (1^+)'' 

{Bhi, -- MdiO'l^l^ ■ • ■ ^O.dXd)}x = (xi,--- ,xa)eR'i = {af^ • • • O-^'^Bhi, -- Mdi^l^ • ■ • iXd)}x& 
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In particular, if we set a = ai = ■ ■ ■ = ad, we recover that 

that is Bhi.- - ,Ha satisfies Property (|2.ip with E = Id and H — Hi + ■ ■ ■ + Hd- 

In Q the existence of OSRF with stationary increments in the stable case 
for any d x d matrix E with positive real parts of the eigenvalues is proved. A 
special class of OSRF is defined through its harmonizable representation. For 
Gaussian models, which is here the case of interest, it reduces to consider an 
integral representation of the form 



for any norm || ■ || on M'^. Such a function / is called a spectral density. In 
order to recover the scaling property (|2.ip , the spectral density / is required to 
satisfy additional specific homogeneity properties (see Section [2] below). In Q, 
such spectral densities are defined by an integral formula. This is a non explicit 
definition in the sense that actual computations require numerical approxima- 
tions. However, these calculations are, in practice, quite difficult to implement. 
Nevertheless, a simpler and explicit formula is furnished in the particular case 
of diagonalizable matrices. 

In this paper, we mainly aim at providing a complete description through 
explicit formulae for the spectral densities in the model defined in Q. We 
focus on a specific case : The Gaussian model. The motivation of this re- 
striction is twofold. On the one hand, it is a reasonable assumption in many 
applications; on the other hand, to improve the model, it is necessary to under- 
stand and classify its geometrical properties which is easier in the Gaussian case. 

Our main results are stated and proved in Section[31 The first ones, Lemma [5Tl 
Lemma l3.21 Lemma 13.31 and Lemma 13.41 

1. reduces the construction of an explicit example of OSGRF for a fixed 
matrix E and an admissible Hurst exponent H (as defined in Section [2]) 
to four particular cases related to specific geometries, 

2. provides an explicit example in each of these four specific cases. 

Thus, we are able to provide an explicit example of OSGRF satisfying Equa- 
tion (jl.ip and then extend the already existing results. Moreover, our second 
result Theorem 13.21 gives a very simple relationship existing between all pos- 
sible spectral densities associated to a given exponent E. This result is not 




where / is a positive valued function defined on M'' 



satisfying 
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formal and can also be turned into an algorithm which generates different fields 
- with different geometries - satisfying Equation (jl.ip for the same matrix of 
anisotropy E. 

These results have important consequences. Firstly, it allows to define the 
studied class of OSGRF from four specific cases. Furthermore we give a complete 
description of the whole class of spectral densities of these fields. Finally, since 
our construction is explicit, the numerical simulations of OSGRF become much 
easier. Thus, our approach provides an explicit definition of an interesting and 
large class of fields for simulations of textures with new geometries. There is 
actually a practical motivation to be able to compare natural/real images (of 
clouds, bones,...) and models with controlled anisotropy. 

In the following pages, we are given d G N \ {0} and E a d x d matrix with 
positive real parts of the eigenvalues. We define 

Xmin{E) - min (Re(A)) . 

\eSp{Ey 

For any a > recall that is defined as follows 

= exp(i;iog(a)) = ^ . 

fe>0 

As usual, i?* denotes the transpose of the matrix E. 

We denote | ■ | the Euclidean norm defined for any x = [xi, ■ • ■ , Xd) € K'' as 



2. Presentation of the model : Operator Scaling Random Fields (OSRF) 

Let us recall some preliminary facts about Operator Scaling Random Fields 
(OSRF) and Operator Scahng Gaussian Random Fields (OSGRF). We refer to 
[6| for all the material of this section. 

Definition 2.1. A scalar-valued random field {X{x)}r^^^d is called operator- 
scaling if there exists a dx d matrix E with positive real parts of the eigenvalues 
and some H > such that 

{X(a^x)}^gRd = {a"X{x)}^^if^d , (2.1) 

(£) 

where — denotes equality of all finite- dimensional marginal distributions. Ma- 
trix E and real number H are respectively called an exponent (of scaling) or an 
anisotropy, and an Hurst index of the field. 

Remark 2.1. In general, the exponent E and the Hurst index H of an OSRF 
are not unique. 
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Thus the usual notion of self-similarity is extended replacing usual scaling, 
(corresponding to the case where E = Id) by a linear scaling involving matrix 
E (see Figure 1 below). It allows to define new classes of random fields with 
new geometry and structure. 




As said in the introduction, when the matrix E is given, the class of OSRF 
with exponent E may be very general. In Q, for any given admissible matrix 
E, the existence of OSRF with stationary increments is proved using a harmon- 
isable representation. 

Recall that according to or [ll[ , given a stochastically continuous Gaus- 
sian field with stationary increments {X{x)}^^g^d, its covariance can be repre- 
sented as 

E{Xix)Xiy)) = / (e'<-^«> - l)(e-'<^'«> - l)dfiiO+ <x,Qy>, 

JR'' 

where Q is a, dx d non-negative definite matrix and fi a non-negative symmetric 
measure fi such that 

jR-i 

for any norm || • || on R'*. Measure fi and matrix Q are proved to be unique. 
Measure /i is called the spectral measure of {X{x)}xfz]^d. In the case where this 
measure is absolutely continuous with respect to Lebesgue measure, the density 
function of /i is called the spectral density of the Gaussian field {X{x)}^^]gd. In 
this case, the Gaussian field {X{x)}^^^d can thus be represented as 

Xix) ^ [ (e'<--«> - l)f/\OdW{0 , (2.2) 
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with 

/ (lA||e||)/m<+oo, (2.3) 

for any norm || • || on R'^. This representation is then called the harmonisable 
representation of the Gaussian field {X{x)}^^g^d. 

To prove the existence of OSRGF with stationary increments for any ad- 
missible matrix E, a quite natural approach is then to use a harmonisable 
representation of the form (|2.2p . In ^] , the following result is proved. We state 
it only in the Gaussian case : 

Theorem 2.1. Let p a continuous function defined on R*^ with non-negative 
values such that for all x G M''\{0}, p{x) ^ 0. Assume that p is -homogeneous 
that is : 

Va > 0, Ve e p(a^'0 = MO ■ 
Then the Gaussian field {Xplx)}^^^.'' defined as follows 

X,{x)= f (e'<-'«>-l)p(0-^-^dl?(e), (2.4) 

exists and is stochastically continuous if and only if H £ (0, A,„in(i?)). Moreover, 
this field has the following properties : 

1. Stationary increments, that is for any h G R'^ 

{Xp{x + h) - Xp(/i)}^gRd ^ = ^ {Xp(a;)}^gRd . 

2. Operator scaling : The scaling relation 112.1]) is satisfied. 

Remark 2.2. Through this new class of Gaussian fields, even if it is a quite 
general model, we do not describe the whole class of OSRGF with stationary 
increments. 

Remark 2.3. If H e (0, Ami„(£^)); /(f) = p{^)-^"-'^''^^'> is proved to be a 
spectral density in the sense that V2. S^) holds. Moreover, the spectral density of 
the Gaussian field {X p{x)} r^i^^d is f . Observe that f is continuous and satisfies 
a specific homogeneity assumption. This homogeneity assumption is necessary 
for the operator scaling property of the Gaussian field {Xp{x)}^^^d whereas 
the continuity assumption ensures that the field {X p{x)} ^^-^d being defined is 
stochastically continuous. 

The main difficulty to overcome is to define suitable spectral densities 
of this new class of Gaussian fields using continuous, i?*-homogeneous functions 
with positive values. In [3l such functions are called (R'', E^) pseudo-norms. 
They can be defined using an integral formula (see Theorem 2.11 of [6]) : 
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Proposition 2.2. Function p defined as 

piO^ / / {l~C08{<x,r^ e>))-dp{9), 

J So Jo ^ 

is continuous with positive values and -homogeneous. Here Sq denotes the 
unit sphere of M'' for a well chosen norm defined from E and /i a finite measure 
on Sq. 

Remark 2.4. Remark that this formula is not the most appropriate for nu- 
merical simulations since we need to approximate an integral. In what follows, 
we will give simpler examples of (R'',£'*) pseudo-norms in the sense that these 
examples lead to exact numerical computations. We then conclude that we give 
explicit examples of (R'^,E^) pseudo-norms (in the numerical sense). 

Remark 2.5. We also refer to P.G.Lemarie (see [16]) whose definition of(R'^,E*) 
pseudo-norms is slightly different (see Remark \3.1\ below). 

Finally, in the special case where matrix E is diagonalizable, an explicit expres- 
sion is given (Corollary 2.12 of [1]) : 

Proposition 2.3. Let E a diagonalizable matrix with positive eigenvalues 

< Ai < • • • < Ad, 

with associated eigenvectors 

Oi, - ■ ■ ,dd, 

and Ci, • • • , Crf > 0. Then for any t < 2Amin(-E') 

p{x)^ (^c,\<x,e,>r/'^ 

is a continuous, E^ -homogeneous function with positive values. 

In this paper, we aim at extending these results and then describing for any 
given admissible matrix E and Hurst index H, all the spectral densities of this 
model of OSGRF with stationary increments in an explicit way. 




3. Definition of explicit spectral densities of the model 

As has already been said in Section [51 the main difficulty is to define appro- 
priate spectral densities of the model. To this end, we note that the class of 
the spectral density used in [6*1 is intimately related to the class of the so-called 
pseudo-norms defined in [16]. We then explicit the link between two (R'^,£') 
pseudo-norms when the matrix E is given. Thereafter using a Jordan reduc- 
tion, for each matrix E with positive real parts of the eigenvalues, we give an 
explicit example of a suitable spectral density of the studied model. Combining 
these two results, we entirely describe in a explicit way the class of spectral 
densities of the Gaussian fields considered in m. 
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3.1. More about pseudo-norms 

Let us first recall some well known facts about pseudo-norms which can 
be found with more details in [l^. This concept is fundamental when defin- 
ing anisotropic functional spaces since using pseudo-norms allows to introduce 
anisotropic topology on W''. Thus, even if the introduction of this concept is 
not necessary to the definition of spectral densities, it is of great importance 
to relate the notion of anisotropic spectral densities to the concept of pseudo- 
norms. This gives us indeed all the tools of "anisotropic functional analysis" 
to study, for example, the sample paths properties of the fields in anisotropic 
spaces (see 0]) and to better understand the inherent topology of these spaces. 

Definition 3.1. A function p defined on R'' is a (W^^E) pseudo-norm if it 
satisfies the three following properties : 

1. p is continuous on , 

2. p is E -homogeneous, i.e. p{a^x) = ap{x) Vx £ M'*, Va > 0, 

3. p is positive on M'* \ {0}. 

Remark 3.1. Our definition of {M.'^,E) pseudo-norm is a slightly modified ver- 
sion of the concept of pseudo-norm on iW'-^A) defined by P.G.Lemarie in fldl l. 
In lid], A denotes a matrix with eigenvalues having a modulus greater than 
one. A pseudo-norm on (R'', A) is a function satisfying properties 1 and 3 of 
the previous definition and the following property : 

p{Ax) — I det{A)\p{x), for any a: m R'^ . 

Further for any dxd matrix A with eigenvalues having modulus greater than one 
and any compactly supported smooth function (j), an example of pseudo-norm on 
(R'', A) is provided by 

=^|det(A)P0(AJx) . 

Remark that if p is a (R*^ , E) pseudo-norm then p(-)^/^'''-"^^ is a pseudo-norm 
on (R'^jA) in the sense of ila] with A = for any given a > 0. The properties 
satisfied by (M.'^,E) pseudo-norms are very similar to those of pseudo-norms on 
(rKa), as proved in m 

Moreover, the example of pseudo-norm on (R'', A) 
given in [16] can be adapted to our case. Indeed for any compactly supported 
smooth function (j) 

p+oc 

P4>{x) = I (j){a^^x)Aa , 

JQ 

IS a (R'', E) pseudo -norm. This formula also leads to numerical approximations 
and thus is a non explicit one. 

The term of pseudo-norm is justified by the following proposition which is proved 
for instance in [l^ or [1] : 

Proposition 3.1. Let p a [W^, E) pseud o-norm. There exists C > such that 
p{x + y)<C{p{x)+p{y)), Va;, yGR''. 
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3.2. Relationship between two given pseudo-norms 

The main result of this section is the description of all the {W^,E) pseudo- 
norms for a given matrix E : 

Theorem 3.2. Let pi be a (R'',^;) pseud o-norm. Then p2 is a (M'', E) pseudo- 
norm if and only if there exists a continuous and positive function g defined on 
E'* \ {0} such that 

P2i0^9iPiiO-''OPi{0- (3.1) 

Proof. Let pi and p2 be two (M'*, E) pseudo-norms. Then the function g — ^ 
is continuous, positive on M"* \ {0} and satisfies for all a > 0, 

g(a^O = 9(0- 

In particular, for a fixed ^ and a = it follows that 

P2(0 = = 5(Pi(0"''0pi(0 • 

The converse is straightforward. ■ 



Consider now the special case E = Id. Theorem 13.21 implies the following 
corollary 

Corollary 3.1. Let {X {x)} r^^jf^d be a Gaussian field with stationary increments 
admitting a continuous spectral density. Assume that X is self-similar with 
Hurst index H . Then, there exists a continuous function S defined on the unit 
sphere £ R*^, |^| — 1} with positive values such that 

(A- w).* ^ {I {^^) s (li) me}^^^^ . (3^2, 

Proof. By assumption the Gaussian field with stationary increments {X{x)}^^]^d 
admits a continuous spectral density denoted /. Then 

^ f (J<x,i> 



X{x) 4 / (e^<-.«> _ l)f^/\i)dW{0 



Since X is self-similar with Hurst index iJ, 

{X{ax)}^(zKd = {a" X{x)}^^^d . 

By assumption, 

X{ax)^ f (e'<'^"'«> - l)/i/2(0diy(0 



Set now C = in the harmonizable representation of X and deduce that 
X{ax) = a-'^l'^ f (e^<-.C> - l)f'^\a-\)dW{0 . 
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Wc now identify the two spectral densities of the two Gaussian fields {X{ax)}xtr]Sid , 
{a^ X{x)}x(z]id which are equal in law . It implies that 



that is p{C) = f(^)-^/(H+d/2) jg (s,djd) pseudo-norm. 

We now apply Theorem 13.21 with E — Id. Then any {R'^,Id) pseudo-norm p 
can be written 



a complete description of self-similar generalized Gaussian fields with station- 
ary r-th increments is given. It implies in particular Corollary . The class of 
anisotropic Gaussian field defined by the representation p.2[) has been widely 
studied (see Recently in [l4], Istas has defined an estimator of S using 

shifted generalized quadratic variations. Let us emphasize that if an anisotropy 
E may be known, using Theorem 13. 21 and a fixed (R'^,i?*) pseudo-norm pi (see 
Section 13.31 below*), one can probably define in a similar way an estimator of 
function g defined in (j3.ip . 

In next section, we now define explicit examples of (R'',i?*) pseudo- norms. 

3.3. Explicit construction o/(R'^,-E*) pseudo-norms 

The result of this section is based on the real Jordan decomposition of any 
d X d matrix E. 

Proposition 3.3. Any d x d matrix E can be written, using the real Jordan 
reduction as 



a 



(3.3) 





where (mi, 7712) £ N x N \ {(0, 0)}, with 
1. For all £1 e {1, • • • , mi}, 



1 



0\ 



Ei^ = \iijd or Ei^ 



1 



Vo 



where A^^ S M, 
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2. For all £2 £ {I,-- - ,m2}, 

/Ae, 



E, 



mi +^2 



with Af,, — 



or E, 



nil +^2 



V 



Af, 



J2 = 



1 




V 

^ where {ai^,l3i^) G 



\ 

h 



As a consequence of the real Jordan decomposition, we state the following propo- 
sition : 



Proposition 3.4. The notations are those of Proposition [3^ For any I, de- 
note df: the size of the matrix E^. Assume that for each 1 < f < wi + m2, 
we are given a (W^^E^) pseudo-norm ti. Define the function (p for any ^ = 

(a,-- - ,Wm2)enrJi""'K'*' as 

1 /2 

Then, the function p defined for any G as 

is a (R'*,£'*) pseudo-norm. Further f = p^(2ff+Tr(E)) ^ suitable spectral 
density of an operator scaling Gaussian random field with stationary increments. 



Proof. Let F = p-^EP. Then for any C - (Ci, • ■ ' , Cmi+m2) e Ilfci 



"il+ni2 • 



-1 (a^ki) + ■ ■ • + T-^i+™2(«^"^+'"<™i+™2) 



- («Vf(Ci) 
= a(p(C) 



mi4-ni2 



(Cmi+ni2)) 



It follows that 

The conclusion is then straightforward. ■ 

Let us illustrate Proposition 13.41 through an example : 



Example 3.5. Set 



E 



2 1 
1 



Note that E is a diagonalizahle matrix since it has two different eigenvalues. 
One has E = PDP-^ with 



D 



2 
\j ' 



P 



1 -1 
1 



11 



A (R'*,_D) pseudo-norm can be defined as 

PD(e) = iCir/' + ie2|. 

Hence Proposition \S.4\ allows to give an explicit expression of a {W^, E^) pseudo- 
norm : 

Ps(0 = pi^(n)-ICir/' + 16-61. 

Remark that in this case, Corollary 2.12 of l^] exactly yields the same result 
since it gives an explicit example of (R'^ji?*) pseudo-norm in the special where 
matrix E is diagonalizable. 

Thus, it is sufficient to define an explicit pseudo-norm for tlie four following 
matrices. 



1. EiiX) 



2. E2{X) = 



3. Esia,^) 



4. E^{a,/3) = 







c 








/A 1 






'■. 1 




V 






vo 


:) ' 



with A 



a (3 
-P a 



(a, 13) e m; X 



/A h 



h 

A) 



with A = 



a 



-P a 



(a,/3) e m; X 



We emphasize that Proposition 13.41 above has two important consequences : 



• The first consequence is that, Lemmas 13. li 13.21 13. 3[ 13^ Proposition 13.41 
and Theorem 13.21 give a complete description of the spectral densities and 
then of the class of Gaussian fields introduced in 6] . 

• Moreover, it implies that all the Gaussian fields belonging to the class 
being studied can be generated from four generic cases corresponding 
to four specific geometries. 

In the four following lemmas, we define a {W^,E) pseudo-norm in each 
generic case. Recall that we denote | • | the Euclidean norm on R''. 



Let us first consider the case E — Ei{\) for some A G R*_: 
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Lemma 3.1. The function pi, defined for ^ gR'^ by 

= i^r^" , (3.5) 

is a (M'^,-Bi(A)) pseudo-norm. 

Proof. The conclusion is straightforward. 
■ 

We now consider the case E = E2{X) for some A G M!J_: 

Lemma 3.2. Let us define the functions Ti and $i for any i e {!, - ■ ■ ,d} as 
follows 

• If 1 = 1, for any £ = (6. ■■■ n{0 = = 161- 

• Ifi>2, for any ^ = (^i, ■■■ ,£d)eM.'' 

'\£^\ ei=6 = ---=C^-l=0, 

. otherwise . 



rm = 

and 



MO = \rim + --- + \nm ■ 

Then, the function p2 defined for ^ G R'' by 

P2{0 = MO'^^ , (3.6) 

is a (R'^jEliX)) pseudo-norm. 

Proof. Let us first prove that the function p2 is well-defined and positive on 
R'^ \ {0}. It is clear that <^d > 0. Further <I>rf(^) = if and only if for any 
i £ {!,■ ■ ■ ,d}, Ti{^) = 0. By induction and by definition of n, it implies that 

iai = --- = ie<ii = o, 

that is ^ = 0. Therefore, p2 is well-defined and positive on R'^ \ {0}. 

To show that p2 is continuous, the only point to verify is that for all 1 < i < d, 

Ti is continuous. To this end, observe that for all ^ e R"^ 

By definition of the exponential of a matrix, one has 



where N = X-^E*{X) - Id. Since N = X'^ 



fe=i 

/O 0\ 
1 ■•. 



one has 



Vo 1 Oj 

(-l)-^log'-^(<J>._i(0), 



-i(0 ($.-i(0 ^), = ^^ + L. (i-ty.xi-i ■ 
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Further one has by induction on i G {2, • • • , d}, that, for ah 1 < ^ < i — 1 

lim(^,log^-^(<i>,_i(a)) = 0. 
Then the continuity of follows. 

We now verify that p2 satisfies the homogeneity condition. It can be done by 
induction on i, showing that, for all 1 < i < d, one has 

r,(a-^^(^)0 = and p.(a-^^(^)0 - 

Indeed, assume that the result holds for i — 1, then for any a > and any ^ such 
that ■ are not both equal to (the other case being trivial), one has 

We have then proved the homogeneity property of function p2 ■ 
■ 

We now consider the case E = E^{a, (3) for some (a, /3) G M^j. x M: 
Lemma 3.3. The function defined for ^ G M'' by 

P3(o = i^r^" , (3.7) 

is a (R'^, Ei^{a, (3)) pseudo-norm. 

Remark 3.2. The function p^ defined by 7] j is an isotropic {W\El{a,l3)) 
pseudo-norm, that is invariant by any isometry T of . Up to a multiplicative 
constant, it is the unique one. Indeed, let p be another isotropic {W^ , E*^{a, (5)) 
pseudo-norm. Observe that for any a > 0, a^a'-"'^-' = a°'T with 



f cos(/3 log(a)) — sin(/3 log(a)) 
I sin(/31og(a)) cos(/3 log(a)) 



fR 
T = ■ . ^ I where R 

VO R 

Remark that T is an isometry. By assumptions on p^ 

P3{a"0 = P^iT-^a^k) = P^ia^'k) = apsiO Va > . 

Then p^ is a (R'*, aid) pseudo-norm. Now, consider the function g defined for 
any ^ G R'^\{0} by g — |CI ""^^"P- Since p^ and \ ■ \^/'^ are two isotropic (W^, aid) 
pseudo-norms, g is isotropic and we have, for all ^ G M.'^ and a > 0, 

9{a"i) - g{0 ■ 
that is setting b — a"' , for all b > 

gm = 9(0 ■ (3.8) 

Hence g is constant on M.'^ \ {0}. 
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Remark 3.3. Using Theorem \3.Sl with g non trivial (i.e. non constant on 
the isotropic unit ball) and pi — \ ■ 1^^", we are able to define non-isotropic 
{R'^, El{a, P)) pseudo-norms. 



Proof. Observe that for any a > 
/i? 0' 



cos(/31og(a)) — siii(/31og(a)) 
sin(/3 log(a)) cos(/3 log(a)) 



with R = 

\0 Rj 

Since the Euchdean norm | ■ | is invariant by any isometry, and in particular by 

/R 0' 



T = 



Vo 



Rj 



one has for any ^ G M'' 

|a"rC| = |a"ei=a"iei • 
The conclusion is then straightforward. 
■ 

We now consider the case E = £^4(0:, (3) for some (a, /3) G x M 
Lemma 3.4. Let us define the functions Ti and for any 1 < i < d as , 

• If 1 = 1, for all e = (fi, ■■■ $1(0 = ri(e) = (|al' + 

• If I > 2, for allC=iii,--- ,U)eM.'' 

+ 16.1')^ «/a = 6 = • ■ • = 6.-2 = 



and 



$.-1(0 



otherwise 



MO = \rim + --- + \nm ■ 

Then, the function p4 defined for all ^ G K'^ by 

PiiO = , (3.9) 

is a (R*^, i?4(a, pseudo-norm. 

Proof. The proof is similar to this of Lemma l3.2l Indeed, set for all 1 < i < d/2 

r.(C) = (|6.-iP + l6.P)", 
and observe that for all ^ G M'^ 

P4{0 = P2{ri{0, ■ ■ ■ ,^d/2(0) > 
where p2 is defined by (|3.6p with A — a. ■ 
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3.4- Two dimensional examples 

We now focus on the two dimensional case. Up to a change of basis, E is a 
matrix of the form : 

1. £;i(Ai,A2) = (^^^ with (Ai,A2) e (m;)2. 

2. E2{X) = with A e r;. 

3. E3{a,P) = (^^^ ^ with {a, 13) £ M^. 

Let us remark that in dimension 2 there is not four generic cases but three since 

/A h 0\ 



the matrix E cannot be equivalent to 



for some (a, (3) G 



\0 A) 



with A = 



a (3 
-(3 a 



We now give an explicit example in each of the case above using the results 
of Lemma 13.11 Lemma 13.21 and Lemma 13.31 

1. If £: Si(Ai,A2) for some (Ai,A2) G {M.*+f, the fimction = 

+ 161^''^')^/^ is a {M?,E\) pseudo-norm. 

2. If = £;2(A) for some A e the function ^2(6, 6) = (161 + 16 - 
4^1n|^i|)i/^ is a (M^,£'2) pseudo-norm. 

3. li E ^ E3{a,l3) for some {a, (3) e the function /9(6,6) = 16^^" is a 
(R'^, £'3) pseudo-norm. More interesting is the fact that the function 

MiCOs(^Mn(r(6)) - 6sin(/3/aln(r(6)| 

P3 = 



with 

is also a {M.'^,El) pseudo-norm. 

Combining these results with Proposition [331 yields us to an explicit example of 
{M.'^,E*) pseudo- norms for any 2x2 matrix E whose eigenvalues have positive 
real parts (see Figure 2 below) . Thereafter Theorem 13.21 brings us a complete 
description of the whole class of (R^,i?*) pseudo- norms for any matrix E and 
thus for spectral densities of the class of OSGRF defined in Q. 



In Figure 2 just below we represented the pseudo-norms pi,p2,P3 for some 



values of the matrix E. Remark that the cases E= il ^] andE ' ^ ^ 



belong to the same generic case (the first one). 
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E=[\ '\ E= ] -}\ E= \ ,M E='' ' 



1/ \1 I ) \Q 1/2) \0 1 



Figure 2: 

Four pseudo-norms corresponding to the three generic two-dimensional cases. 
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